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Two electron spins in quantum dots coupled through coherent tunneling are generally acknowl- 
edged to approximately obey Heisenberg isotropic exchange. This has not been established for two 
holes. Here we obtain the spectrum of two holes and of two electrons in two vertically stacked 
self- assembled InAs quantum dots using optical spectroscopy as a function of electric and magnetic 
fields. We find the exchange to be dominantly isotropic for both systems, with observed anticrossing 
singularities that can only be explained by asymmetric mixing. The latter arises from spin-orbit 
and Zeeman interactions combined with spatial asymmetries. Such asymmetric contributions to the 
isotropic Hamiltonian for electrons are of the order of a few percent while those for holes are an 
order of magnitude larger. 

PACS numbers: 71.70.Ej, 71.70.Gm, 78.55.Cr, 78.67.Hc 



The exchange interaction between the spins in two 
quantum dots (QDs) leads to entanglement and to the 
opportunity of quantum information processing |l| . Elec- 
tron spin qubits have been studied in great detail, but 
recently hole spins have received special attention 0, Q 
because of their reduced hyperfine interaction with the 
nuclear spin reservoir Now, for inter- QD entangle- 

ment, the nature of the exchange interaction is of central 
importance, both for quantum gates and for decoherence 
of two-qubit states. 

The symmetry of this interaction is often assumed to 
be the Heisenberg isotropic exchange, J<j\ • CT2- As a re- 
sult, the eigenstates of two spins in two quantum dots 
form singlet and triplet spin states separated by the ex- 
change energy, J. This model is very important, concep- 
tually [l| and for interpretations of complex experimen- 
tal spectra jgj], and has been used widely in magnetism, 
quantum computing, molecular and quantum dot struc- 
tures. It has been shown that for two bound electrons 
the isotropic Heisenberg interaction captures almost all 
of the physics and it requires only small asymmetric ex- 
change terms that couple singlets and triplets (9l-[l4j. 

These small singlet-triplet mixings are important in 
spin relaxation in bulk GaAs |15j. In QDs, they require 
solutions to remove their effects on the fidelity of two- 
qubit gates [Irl[l7|. 



In contrast to electron spins, hole spins are extremely 
anisotropic, for example, in their ^-factor [18j and in 
their hyperfine interaction [19|- This anisotropy arises 
from the strong spin-orbit character of the valence band. 
Counterintuitively, we find that isotropic hole-hole ex- 
change between two self-assembled InAs QDs is in fact a 
good starting point [lo|. This surprising result was pre- 
dicted by theory, which demonstrated that in the absence 
of heavy hole -light hole mixing, exchange between two 



> 



O) 

CD 
C 
LU 



B = 




B > 0T / 

/ T. 






2V2ti/ 


^^r~s(i,i) 


t S 


/ T+ 

/S(2,0) (b) 



Electric Field, F Electric Field, F 



Magnetic field, B 



FIG. 1. (a) Hund-Mulliken model for two-spin configuration 
at zero B. (b) F dependence at fixed B for a symmetric case, 
Eq. ©. (c) B depend ence at fixed F (marked by arrows with 
Jin(b)). 



hole spins in two QDs should be isotropic, just like two 
electron spins [2l|. Now we find that substantial asym- 
metric contributions arising from spin-orbit interactions 
and inhomogeneous Zeeman interactions are necessary to 
understand the measured 2h energy levels. 

The asymmetric contributions are often overlooked in 
spectroscopy due to their small perturbations to QD en- 
ergy levels, where the energies are dominated by the 
isotropic exchange. In this paper we present an optical 
study to identify and measure asymmetric contributions 
to the interaction between two holes (2h) and also be- 
tween two electrons (2e) in two tunnel-coupled QDs. We 
use individual pairs of vertically stacked self-assembled 
InAs/GaAs QDs separated by a thin tunnel barrier with 
a thickness d. Two types of samples were developed using 
a Schottky diode grown by molecular beam epitaxy, one 
for 2h [3] and one for 2e (22[. The width and height of the 
tunnel barriers were selected in a way that allowed us to 
achieve similar values of singlet-triplet splitting for both 
cases (J « 100 /ieV) Q. Electric (F) and magnetic (B) 
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fields were applied longitudinally along the stacking z- 
axis in the Faraday geometry and are used to balance the 
dominant isotropic exchange and to bring the energy lev- 
els to resonance displaying their mixings. Measurements 
were done at 5 K using photoluminescence with a spectral 
resolution limited by the spectrometer of ~ 15 /ieV, and 
laser transmission spectroscopy with a resolution deter- 
mined by the linewidth (~ 6 fieV for holes and ~ 3 fieV 
for electrons). 

The Heisenberg exchange can be treated within the 
Hund-Mulliken model 0, HH. The natural spin state 
basis is three triplets: To = (t, ^) To , T + = (t?t)T+> 
T_ = (|, l) T _ and three singlets: S^,o) = CNoO)s> 
5(1,1) = (t4) 5 , and 5(0,2) = (0,tl) 5 . The individ- 
ual spin projections are either the electron spin ±1/2 
or hole pseudo-spin ±1/2. The singlet states are cou- 
pled together and shifted in energy by spin conserving 
tunneling (t) between the two QDs, but because of spin 
blocking, the triplet states are not affected. The Hamil- 
tonian within the singlet basis is given by: 



(tl,0) 5 (ti) s (0,tl)s 

'Fd+U -y/2t 

-y/2t -V2t 

-y/2t -Fd+U, 



(1) 



The potential (U) is the Coulomb energy required to 
move the two charges from separate QDs to the same 
QD. The relative energy between QDs separated by d is 
controlled via F. In Fig. [1Ja), the resulting levels of the 
singlets show anticrossings. The lowest energy singlet 
state is S = aS( 2j o) ± ^(i,i) + c ^(o,2)- In Fig. QJb) and 
hereafter we focus on one of the anticrossings, for which 
c = 0. The isotropic exchange interaction J is defined as 
the splitting between To and the lowest singlet with the 
spin Hamiltonian Jcri • (T2, with J = J(T, U, t). 

To fully probe and engineer the spin states of the two 
QDs we also need, in addition to T, a magnetic field. 
The simplest Hamiltonian consists of isotropic exchange 
and an average Zeeman interaction, 



H = Ja ± - cr 2 + P ext -(0-1+0-2) 



(2) 



The Zeeman term (3 ext = (gn ± #22)mb^/2 splits the T + 
and T_ lines, but preserves the spin states. In Fig. QJc) 
the energies are calculated as a function of B with F held 
constant. We find experimentally that this symmetric 
approximation is fairly good for two electrons but less so 
for two holes. 

Calculated transitions from the 2h levels to the charged 
exciton levels (X 2+ ) are shown within the symmetric ap- 
proximation in Fig. [2fa). Here we show the dependence 
on B, holding F constant. The exciton states have been 
described previously @, [HI . The corresponding transi- 
tion spectrum is shown in Fig. [2^b). At zero B the tran- 
sition spectrum shows a doublet with a splitting given by 
the exchange J, which then separates into four transition 
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FIG. 2. (a) Calculated B dependence for the excited and 
ground states for the 2h in symmetric case. Transitions are 
labeled by final state, (b) calculation of the transition spec- 
trum for the B dependence of X 2+ using the parameters: 
9e = -0.46, g h ,! = g h , 2 = 1.1, t = -115/xeV, U = lOmeV. 
The diamagnetic shift is not taken into account. 



lines with increasing B. Transitions involving £ are non- 
degenerate at finite £?, but those involving To, T± remain 
two-fold degenerate - in the absence of asymmetry. 

In contrast to the symmetric calculation of Fig. [2fb), 
the measured 2h spectrum appears much more complex 
as shown in Fig. [3fa). In addition to substantial line 
splittings observed in the triplet optical transitions that 
grow with increasing 5, there are now anticrossings with 
magnitudes of 26 fieV at B = 1.5 T. These anticrossings 
cannot be captured by a symmetric hamiltonian and re- 
quire off-diagonal matrix elements that describe coupling 
between S and T. We will call them "spin-flip" anticross- 
ings. We have also measured the 2h optical spectrum as 
a function of F at fixed B. The complex pattern of an- 
ticrossings observed in the data of Fig. [3jd) arises from 
spin-conserving anticrossings of both the 2h and the ex- 
citon states as discussed previously jlH combined with 
the smaller spin-flip anticrossings. All of these features 
can be explained with natural extensions of the symmet- 
ric Hamiltonian Eq. (j2j) as we discuss below. As a result 
we obtain the fitted optical transition spectrum shown 
in Figs. [3fb) and[3fe), and calculate the corresponding 
2h energy levels in Figs. [3fc) and [3^f). A comparison 
of these results and the symmetric calculations of Fig. Q] 
is additional evidence that asymmetric contributions are 
important. 

The 2e optical spectrum found in a manner similar to 
that of Fig. [3fa) gives a much smaller spin-orbit contri- 
bution than for the 2h case. In fact it was necessary to 
use the higher resolution of laser spectroscopy to resolve 
the additional structure [22| . The magnitude (4 /ieV) of 
the spin-flip anticrossing was 4% that of the symmetric 
exchange (100 fieV) and 1% that of the spin-conserving 
tunneling anticrossing. In addition, the splitting in the 
triplet optical transition at high fields (4 //eV at 4 T) due 
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FIG. 3. (a) measured and (b) calculated transition spectra 
of B dependence for the X 2+ . gh,i = 1.85, gn,2 — 0.78, 
9e = -0.46, g h , 12 = 0.14, t = -115 /leV, 712 = 24/leV, A 7 - 
5 /xeV. (a) and (b) Three signatures of AC are marked by 
dashed horizontal arrows. A diamagnetic shift of 11.6 /xeV/T 2 
is subtracted in (a), (c) Calculated 2h energy levels with 
ground state S/T- AC of 26/xeV. (d) measured and (e) calcu- 
lated transition spectra of F dependence for B — 9 T. (small 
black arrows on the right of panel (d) mark a transition spec- 
tra of X + .) (f) Ground state energy levels as a function of 
bias. Circles mark the spin-flip ACs. Vertical lines mark the 
position of ACs in transition spectra of Fig. [3{d) and [3{e). 



to ^-factor inhomogeneity was small. 

In the data modeling shown in Fig. [3] we have gener- 
alized the symmetric spin Hamiltonian of Eq. (J2j) to in- 
clude spin coupling terms as follows. The spin interaction 
of holes (or electrons) again is included with interaction 
terms, f3\ ot • <j\ + ft 1 ^ 1 • <J 2 acting on each spin. The fields 
now include not only the external B but also the inter- 
nal relativistic magnetic fields arising from the spin-orbit 
interaction (3\ ot = (3^ xt + (3f°. We could also include the 
Overhauser fields due to the hyperfine interaction with 
nuclear spins, in a similar way to Ref. [14[, but this is ex- 
pected to be small here. After integrating over the orbital 
degrees of freedom and including the isotropic exchange 
we obtain the general spin Hamiltonian for two spins in 
two QDs: 



H = JCT! • <7 2 +f3 - (CTI + <7 2 ) + • (CTI - <T 2 ) + P" ' Ol X <7 2 ) 

(3) 

In general each of the interaction terms can depend 



on both external and internal fields. In our case, the 
strong axis of the QDs and also the external B are along 
the z-axis, B = Bz. The QDs are also stacked along 
the z-axis, and we take any lateral asymmetry (e.g. an 
offset between the centers of the QDs) to lie along the x- 
axis. As a result the internal field (3 so = jx. With these 
Zeeman and spin-orbit fields the parameters in Eq. (j3j) 
can be written: 

(3 = l/2fj jB B? l gz + Zjx 
f3 f = l/2fi B B(ag 12 + bAg)z 
P" = O712 + bAj)x. 

The Zeeman interaction with the external B is de- 
termined by the matrix elements of the coordinate- 
dependent ^-factor g(r) over the two QDs: £g = (gu + 
022), A# = (011-022) and 012 = <l|^(r)|2>. Likewise, the 
spin-orbit field leads to mixing terms given by matrix el- 
ements of 7. The linear superposition aS^ y o) + ^(1,1) 
gives a contribution of 712 and A7, which is determined 
by measuring it's dependence on F. 

We can also write the Hamiltonian of Eq. (j3j) in a ma- 
trix form within the singlet /triplet spin basis that shows 
explicitly how the states are mixed, and permits conve- 
nient fitting to the data. 
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(4) 

Using Eq. ((4]) we were able to get good fits to the data 
as shown in Fig. 02 We obtain the coupling parameters 
pg,Ag,g 12 ) = (2.63,1.07,0.14), and (£7^7,712) = 
(0,5,24) /ieV. 

The diagonal terms correspond to the symmetric 
Hamiltonian of Eq. (|2j). The new off-diagonal terms lead 
to the observed energy splittings and anticrossings in the 
spectra. The phenomenological Zeeman and spin-orbit 
parameters, g^ and 7^, have useful physical interpreta- 
tions. First, the term with Ag is manifestly the difference 
in factors between the two QDs, and likely arises from 
the difference in size and indium concentration between 
the two QDs. This parameter is larger for holes than 
for electrons in part because the larger effective mass of 
the hole makes them more localized and more sensitive 
to differences of the QDs. The term with g\ 2 physically 
arises from the difference in ^-factor between the barrier 
and the QDs, which leads to the t and \, spins-tunneling 
at different rates [26]. This terms breaks the spin sym- 
metry and mixes S and To. At sufficiently high fields the 
eigenstates become (t, I) and (^, t), instead of (t, i)s and 

(tl)To 0- 
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These off-diagonal Zeeman terms effectively push the 
S and T energies apart as a function of B as seen in 
Fig. [3^c), arrows on the right. This is measured directly 
in the optical spectrum by the splitting in the triplet 
transitions that grows with B. This splitting is poten- 
tially important because it enables resonant optical con- 
trol and measurement of each of the spin states sepa- 
rately. Recent studies have demonstrated optical quan- 
tum control of coherent superposition of the To and S 

I 1 I li , I L 



for 2e [22|, [28| and 2h [3[ states. New opportunities in 



control and measurement may be possible with further 
increases of the triplet splitting that arises from differ- 
ences in ^-factors. For example, if the S and To states 
are used as an optically controlled qubit, the T_ and T + 
transitions could provide cycling transitions for readout, 
analogous to what has been done with a single electron 
qubit system [29[. 

The result of the spin-orbit field is similar but with an 
important difference. The spin-orbit field acting orthog- 
onal to z has the effect of partially rotating t to ^, and 
vice versa. In particular, the term 712, which couples 
(11,0)5 O (tjt)T_j can be viewed as spin- flip tunneling. 
It is analogous to the spin conserving tunneling term (t), 
which couples (1\U0)s ^ (t?i)s- The 712 term along 
with the A7 term have the effect of mixing the S state 
with the T_ and T + triplet states, and leads to the an- 
ticrossing observed in Fig.[3fc). Quantum control of two 
holes in this system could thus be obtained with elec- 
trostatic gates in a way analogous to Taylor et al. [30| , 
however, instead of hyperfine coupling, the much larger 
spin-orbit interaction could be used. 

The term £7 couples To with the T + and T_ triplets. It 
can be measured as a splitting of the triplet line at zero 
B. Any zero-field splitting of the triplet energies was 
found to be less than our resolution of 15 jaeV for the 2h 
transitions, and so we took this parameter to be zero [3l|. 
We note that in a separate study there is evidence for a 
very small splitting of ~ 8 /aeV in a similar sample from 
coherent measurements in the time domain [3[. 

Microscopically, the f3 so field arises from spin-orbit in- 
teractions combined with spatial asymmetry (e.g. the 
lateral displacement of the two stacked QDs). It does 
not couple the states directly but rather results indirectly 
from coupling with excited states [32]. More generally, 
the magnitude of asymmetric exchange will depend on 
the width of the tunneling barrier and on the magnitude 
of spatial asymmetry (e.g. misalignment of the QDs). 

Here we have found that a symmetric spin Hamilto- 
nian based on the isotropic Heisenberg exchange interac- 
tion can be generalized to treat the 2h as well as the 2e 
spectrum in tunnel-coupled QDs using phenomenological 
off-diagonal Zeeman and spin-orbit parameters. 



[1] 
[2] 

[3] 
[4] 

[5 
[6] 

[8 

[9 
[10; 

[11 
[12 
[13 
[14 
[15 

[16 
[17 
[18 
[19 
[20; 



D. Loss and D. P. DiVincenzo, 
|Phys. Rev. A 57, 120 (1998)] 

K. De Greve, P. L. McMahon, D. Press, T. D. Ladd, 
D. Bisping, C. Schneider, M. Kamp, L. Worschech, 
S. Honing, A. Forchel, and Y. Yamamoto, 
|Nat. Phys. 7, 872 (2011)] 

A. Greilich, S. G. Carter, D. Kim, A. S. Br acker, and 
D. Gammon, |Nat. Photon. 5, 702 (2011)] 

D. Brunner, B. D. Gerardot, P. A. Dalgarno, G. Wst, 
K. Kar rai, N. G. Stoltz, P. M Petroff, and R. J. War- 
burton, |Science 325, 70 (2009) | 

B. Eble, C. Testelin, P. Desfonds, F. Bernardot, A. Baloc- 
chi, T. Amand, A. Miard, A. Lemaitre, X. Marie, and 



M. Chamarro, Phys. Rev. Lett. 102, 146601 (2 009) 
P. Falla hi, S. T. Yilmaz, an d A. Imamoglu, 
|Phys. Rev. Lett. 105, 257402 (2010) | 
E. A. Chekho vich, A. B. Krysa, M. S. Skolnick, an d A. I. 
Tartakovskii, |Phys. Rev. Lett. 106, 027402 (2011)1 
M. Scheibner, M. F. Doty, I. V. Ponomarev, A. S. 
Bracker, E. A. Stinaff, V. L. Korenev, T. L. Reinecke, 
and D. Gammo n, |Phys. Rev. B 75, 245318 (2007)] 
K. V. Kavokin, [Phys. Rev. B 64~075 305 (200ljp 
S. C. Badescu, Y. B. Lyanda-Geller, and T. L. Reinecke, 
[Phys. Rev. B 72, 161304 (2005) | 
S. 



Chutia, 



M. Friesen, 



and R. Joynt, 
and J. Fabian, 



Phys. Rev. B 73, 241304 (2006) 

F. Baruffa, P. Stano, 
|Phys. Rev. Lett. 104, 126401 (2010)] 
M. P. Nowak an d B. 
|Phys. Rev. B 82, 165316 (2010) | 
D. Stepanenko, M. Rudner, B. I. Halperin, 
|Phys. Rev. B 85, 075416 (2012) | 

R. I. Dzhioev, K. V. Kavokin, V. L. Korenev, M 
Lazarev, B. Y. Meltser, M. N. Stepanova, B. 
Zakharchenya, D. Gammon, 
|Phys. Rev. B 66, 245204 (2002) | 

G. Burkard, D. Loss, and 
|Phys. Rev. B 59, 2070 (1999)] 
G. Burkard, G. Seelig, 



Szafran, 



and D. Loss, 



V. 
P. 



and D. S. Katzer, 
D. P. DiVincenzo, 



and D. Loss, 



Vogl, 



Loss, 



Phys. Rev. B 62, 2581 (2000) 

T. Andlauer and P. 

|Phys. Rev. B 79, 045307 (2009)) 

J. Fischer and D. 

|Phys. Rev. Lett. 105, 266603 (2010) | 
As a counter example, electron-hole exchange in the case 
of the neutral exciton is nearly Ising-like (Ja ez <Jhz) with 
two nearly degenerate doublets separated by the ex- 
change splitting. In that case, a small asymmetric e-h 
exchange 8(0^0 y + cfy(T%) mixes the two bright excitons. 
Another example is the case of two holes in a single QD - 
one in the s-shell and the other in the p-shell. In contrast 
to two electrons in a single QD, the energy structure de- 
viates strongly from isotropic, and is apparently neither 
Heisenberg- like nor Ising-like [33]. 
[21] K. V. Kavokin, jPhys. Rev. B 69, 075302 (2004) | 
[22] D. Kim, S. G. Carter, A. Greilich, A S. Bracker, and 

D. Gammon, Nat. Phys. 7, 223 (2011)] 
[23] For the hole sample we have used a Be-doped GaAs 
buffer, 25 nm GaAs spacer, 2.8 nm InAs bottom dot, 6 nm 
GaAs tunnel barrier, 3.2 nm InAs top dot, followed by 
280 nm GaAs capping layer. The electron sample was 



5 



similar but with a 9nm barrier composed of 3/3/3 nm 
of GaAs/Ga.7A1.3As/GaAs 22]. Both samples were cov- 
ered with a 120 nm Al shadow mask perforated with one 
micron apertures to enable the measurement of single 
pairs of dots. 

[24] W. G. van der Wiel, S. De Franceschi, J. M. Elzer- 
man, T. Fujisawa, S. Tarucha, and L. P. Kouwenhoven, 
|Rev. Mod. Phys. 75, 1 (2002) | 

[25] M. F. Doty, M. Scheibner, A. S. Bracker, I. V. 
Ponomarev, T. L. Reinecke, and D. Gammon, 
|Phys. Rev. B 78, 115316 (2008)| 

[26] M. F. Doty, M. Scheibner, I. V. Ponomarev, E. A. Sti- 
naff, A. S. Bracker, V. L. Korenev, T. L. Reinecke, and 
D. Gammon, Phys. Rev. Lett. 97, 197202 (2006) | 

[27] For holes, the difference in ^-factor of the InAs QDs 
and the GaAs tunnel barrier is substantial and g±2 — 
0.14 [34]. For electrons it was found that this differ- 
ence was negligible with a GaAs barrier, however, it has 
been found recently that with a GaAs/AlGaAs/GaAs 
barrier, as used here, this term can also be significant 
(gi2 = 0.3) [1^]. The symmetric Hamiltonian by itself 
provides a good approximation to the 2e case. 

[28] K. M. Weiss, J. M. Elzerman, Y. L. Delley, 
J. Miguel- Sanchez, and A. Imamoglu, (2012), 



|arXiv:1112.4710v2 [cond-mat.mes-halTJI 

[29] D. Kim, S. E. Economou, S. C. Badescu, M. Scheib- 
ner, A. S. Bracker, M. Bashkansky, T. L. Reinecke, and 
D. Gammon, [Phys. Rev. Lett. 101, 236804 (2008)| 

[30] J. M. Taylor, J. R. Petta, A. C. Johnson, 
A. Yacoby, C. M. Marcus, and M. D. Lukin, 
|Phys. R ev. B 76, 035315 (2007) | 

[31] We have set all 7™ terms to zero at B ~ in Eq. J4}. 
A7 ^> E7 indicates that the local spin-orbit fields /3l°, 
/3|° in the two QDs have opposite directions and similar 
magnitudes. This does not allow for any splitting between 
the To and the T+/T_ states, which must be less than 
our resolution of 15 /xeV. 

[32] M. F. Doty, J. I. Climente, A. Greilich, 
M. Yakes, A. S. Bracker, and D. Gammon, 
|Phys. Rev. B 81, 03 5308 (2010) | 

[33] M. Ediger, G. Bester, B. D. Gerardot, A. Badolato, P. M. 
Petr off, K. Karrai, A. Zunger, an d R. J. Warburton, 
|Phys. Rev. Lett. 9 8, 036808 (2007) 

[34] M. F. Doty, J. I. Climente, M. Korkusinski, M. Scheib- 
ner, A. S. Bracker, P. Hawrylak, and D. Gammon, 
|Phys. Rev. Lett. 102, 047401 (2009)] 

[35] W. Liu, S. Sanwlani, R. Hazbun, J. Kolodzey, 
A. S. Bracker, D. Gammon, and M. F. Doty, 
|Phys. Rev. B 84, 121304 (2011)| 



